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We derive exact matrix integral representations for different sums over partitions. The character- 
istic feature of all obtained matrix models is the presence of logarithmic (or, vice versa, exponential) 
terms in the potential. Our derivation is based on the application of the higher Casimir operators. 

(-») . The Toda lattice integrability of the basic sums over partitions can be easily derived from the matrix 

fl} ' model representation. 
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Introduction 

Random partitions is a popular subject of modern mathematical physics. They appear in such diverse 
areas as 2d Yang-Mills [TJ-[3] and 3d Chern-Simons theory [4j{6] , instantonic calculus of supersymmetric 
gauge theories in different dimensions [THTU] and Hurwitz-Hodge-Gromov-Witten theory [TH - TH] . Of 
course, sums over partitions/representations play an increasing role in string theory, so all aforementioned 
theories can be described by particular string models [31E1IH1HH1HS] • It is the topological nature of the 
considered in all above mentioned examples invariants what lets one to calculate at least some of them 
in the domain of the topological string theory with its powerful topological vertex machinery |17U18| . 
Recently the subject of the partition/representation summation has appeared in the profound AGT 
conjecture [191120] . which connects the partition functions of some supersymmetric gauge theories and 
related string models with the conformal blocks of 2d conformal held theories. 

Sums over partitions are the discrete analogs of the matrix models. This is a different class of models, 
extremely important for modern theory (for a review of recent developments see [21][22| and references 
therein) . Thus it is not unreasonable to ask a question: what are the precise relations between the models 
of two families? An answer to this question is important for the investigation of the models of both types. 
This especially concerns a less developed theory of random partitions. The point is that identification of 
the sums over partitions with matrix integrals allows one to use more elaborated theory of matrix models 
with its powerful Virasoro constraints, well developed semiclassical techniques, dualities between different 
matrix models and rich integrability properties for less developed theory of random partitions. Of course, 
it is well known how to apply large N matrix models techniques to different sums over partitions [101124) , 
but here we are for exact relations. 

Some exact relations between matrix models and random partitions are very well known. Perhaps the 
simplest examples of such connections are character expansions of the celebrated Itzykson-Zuber matrix 
integral (our notations are explained in Section [1} 

[dUle Tr(UAutB) = V^ ^aXa(A)xa(B) 

^—^ dim\ \ ' 



NxN 



\;l(X)<N 



and the unitary matrix model 



/ [dU]exp(f>TrU fc +< fc TrUtM = £ Xx(t)xx(t) (2) 

JNXN \k=0 J X:l(X)<N 

In both examples unitary matrix integrals are equal to the sums with summands made of Schur functions. 

However, in the majority of interesting applications mentioned above more involved sums over parti- 
tions appear. For example in the sums describing 2d YM or double Hurwitz numbers there appear (the 
eigenvalues of) the quadratic Casimir in the exponential. Higher Casimirs also appear in other examples. 
Of course, it is usually simple to switch on the first Casimir, which counts the weight of the partition 
|A| = J^Aj, but for the higher Casimirs the construction of the related matrix models can be rather 
nontrivial. 

Several important examples of the relations between sums over partitions with higher Casimirs and 
matrix integrals are known. Probably the most illustrative example is the noncommutative U(l) gauge 
theory, which is dual to the stationary sector of type A topological string model on CP 1 . Partition 
function of this model is given by the sum of the random partitions, where higher Casimirs correspond 
to the descendants of the Kahler class [8HTT|. An old-standing conjecture [23] states that the partition 
function is given by the Eguchi-Yang matrix integral. This conjectured integral holds a number of 
important properties of the partition function, however it has not much chance to give the correct answer 
in its simplest form (see [TU] for attempts to refine this matrix model representation). Recently it was 
shown that corresponding sum over random partitions is given by a matrix integral of another type [25] . 
In this last matrix integral the sum over partitions appears naturally as the sum over residues in the 
eigenvalue integral due to special choice of the integration contours and the potential. Another important 
example is the generating function of the simple Hurwitz numbers. Here two different matrix integrals 
are known: one [SSj developing the ideas of [35], and another [37] with usual integration contours but 
non-flat measure. These two matrix models are related through the Fourier-Laplace transform [35] . 



In this paper we derive the matrix integral representations of rather general random partitions models. 
These representations are close in their spirit to integrals discussed in [27] and [4^6], and to some extent 
generalize them. We restrict our attention to the sums of the following forrro 

Z%£"\lP-\ . . . , t« X x , . . . , X,; s) := £ dJdim?XA(* (1) ) ■ • • Xx(t {k) )xx(*i) ■ • • Xx(^l)e^ s " C " 

l(X)<N 

(3) 
where the summation runs over all partitions (Young diagrams) with number of parts smaller or equal 
to N. We explain our notations in Section [TJ here let us just make a few comments. The summand in 
formula ([5]) consist of two parts: the potential and the measure. The potential is made of the eigenvalues 
of the Casimir operators 

C fc -g(A i -i + i)*-(-< + i)* (4) 

As in (JTJ and ((2J) the measure is build of Schur functions that depend either on infinite set of time variables 
t or on N x N matrix X. On the matrix model side these two types of variables are naturally identified 
with two main ways to introduce the coupling constants. Namely, in the simplest case of the integration 
over Hermitian matrices, one can either couple times tk with traces of the matrix powers to construct 
usual Hermitian matrix model or introduce an external matrix as in the Generalized Kontsevich Model. 
Let us emphasize that the dependence on times is more general and universal than the dependence on the 
external matrix. The simple reason for this is that when you know the function dependent of the infinite 
set of times t you can easily restrict this dependence on the iV-dimensional subspace by the Miwa change 
of variables tk = ^Tr X* . To be able to make an inverse change in principle one should consider the matrix 
of the infinite size, and even in this case inverse transformation on the level of the partition function can 
be not so transparent. This obstacle is well known for Generalized Kontsevich Model (GKM) [30! : while 
expressions for the partition functions are known very well in terms of the external matrix, no simple 
expression in terms of the time variables is available (see, however, 1311 for the cubic Kontsevich model). 
Thus we present formulas with dependence on times t when possible, and specify them to Miwa variables 
with a matrix X only when this leads to a significant simplification. Other elements of the sum §3§ are 
the dimensions of symmetric and general linear groups representations labeled by the partition. They 
are particular values of the Schur functions: d\ = x\{tk = ^fe,i) and dim a = xa(1)- 

As we do not know any natural examples of random partitions with the " dynamical" Schur functions 
Xx{t) or xa(X) standing in the denominators, we assume that only n and m can be negative, but k > 
and I > 0. If the sum m + I is nonnegative the restriction on l(X) in the sum ([3]) is excessive, and we 
will freely omit it, otherwise this restriction is required. While, generally speaking, the parameter p 
defined by the constraint 2 — 2p = n + m + k + l can be both integer and half-integer as well as both 
positive and negative, only non-negative integer p fits well into our matrix integral construction. In this 
case p corresponds to the genus of the corresponding target manifold and hereafter we assume it to be 
non- negative integer. 

Of course, not all sums (|3|) are independent. One can find several simple relations between different 
sums of this kind. First of all, time variables can be substituted by Miwa variables tk = rTr X fe for some 
matrix X. Then, one can further specify tk = $k,i or X = 1 in one or several Schur functions to get d\ 
and dim \ respectively. One can further exchange d\ and dim \ by the cost of change of the potential 

Y^j S kCk- 

It is also possible to glue two partition functions (or two Schur functions inside one partition function) 
"along" one of the matrices X with the help of the unitary or complex matrix integrals. Thus to get 
matrix model representations for all partition functions ([3]) with non-negative integer genus p, one needs 



1 Of course, this in not the most general form, which appears in the applications, even if one does not consider sums over 
multiple partitions, which are of primary interest for some applications [8 .10 19 ,20 • In particular, much more involved sum 
describes the full partition function of the CP model |16j . Another important generalization is given by the generating 
functions of the generalized Hurwitz numbers, in which one exponentiate not only Casimirs Cfc, but their polynomial 
combinations, namely profound cut-and-joint operators. This type of sums is a direct discrete analog of ordinary matrix 
models with multi-trace potentials. Let us also mention here very interesting and important /3-deformations 19 , 20 29 
and g-deformations [5.6.17 22 25^. We hope to return to the matrix models for those modifications in the subsequent 
publications. 
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Figure 1: The propagator P/v(Xi, X 2 ; s) 





Figure 2: The vertices W^v(Xi,X 2 ,X 3 ) and VAr(Xi,X 2 ,X 3 ). 

to know the matrix models for the finite number of basic functions corresponding to simple topologies. 
To be more specific, one needs the matrix models for vertices of two types: with d\ 
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in the denominators and the matrix model for the propagator with all Casimirs: 



(5) 



(6) 



Pjv(Xi,X 2 ;s) := Zffi' (Xi,X 2 ', s) = J^ xa(Xi)xa(X 2 ) exp^] s l C l 

1(\)<N i=l 

Vertices dU and ([6]) are given by specific complex 

V N (t,A,B)= f [dZ]exp(-TrZZ t + ^i fc Tr(ZAZ t B) fc j 



and unitary 



W N (t, A, B) = / [dU] exp ( V i fe Tr (UAU t B) fc 
•™ \fe=i 



(7) 



(8) 



(9) 



matrix integrals respectively, while propagator (|7|) is of primary interest for us. Indeed, to construct 
all partition functions ([3]) it is enough to know P/v(Xi, X 2 ; s), while a dependence on times t can be 
restored by gluing with a function P/v(£,Xi j2 ;0) = exp^ifcTrX fc . However, our derivation shows that 
the propagator P^(t,i; s) with the Schur measure [13] is more symmetric and two matrix integrals, that 
is one in P/v(i, X; s) and another which glue it with the function P/v(t, X; 0), unify into one integral over 
normal matrices in the very nice way. Different specifications of this function play the most important 
role for the whole story of random partitions and frequently appear in applications. 



We claim that the most natural language for construction of the matrix model representations of 
([3]) and, in particular, of (JT)) is the language of the Casimir operators Ck- We describe three different 
representations of such operators and derive explicit expressions for all operators Ck- Being exponentiated 
these operators lead to the "matrix integral-valued" differential operators, which we use. 

As an example of our approach we consider propagator, in which only the second Casimir Ci = 
~Y^i A-t(A.j — 2i + 1) is inserted. There exists a huge class of applications where this specification plays the 
main role 2-6,26,28 . To work it out we take only s\ and s-i in ([7]) to be nonzero - we denote them by 
q and | respectively. Then for the propagator ([7]) we get the Hermitian matrix integral with a non-flat 
measure 



P N (t,e" 



f [dfiCY)] exp [ -Tr *Y - ^-Tr Y 2 
h \9 2 9 



N 
~2 



TrY + ^£ fe Tr< 



A:Y 



fc=l 



(10) 



The proportionality constants for this and for the subsequent matrix integrals do not depend on times 
tk and can be obtained from the obvious equality P/v(0, •) = 1. 

We managed to show that for two sets of times propagator is given by the following normal matrix 
integral 



P N (t,t) 



I — [dZ L - cx p (-^ Tr lo § 2 zzt + E (^ Tr z" + ^- Tr ztfe ) 

U (detZZt) A ' + 2-i \ 2g ^ 



(11) 



Then we turn on all coupling constants s& and put them to be the Miwa variables Sk — j;Tr Y fe . In 
this case instead of (|10[) we get a complex matrix integral 



P N {t,e*;Y) - / [dZjexp^TrZZ+Y + ^ZtZH-*)) 



with the potential 



N 
#(A) = --TrA + ^i fc TY ( ' 

fe=i 



frA 



E 



(-iy b. 



i+j 



i;j=0,i+j>0 



2(i + j) i\j\ 



TrA'TrA J 



(12) 



(13) 



where Bk are Bernoulli numbers. We want to stress here that contrary to (fTU)) . (fTTT) and (fT4l) this matrix 
model is not immediately reducible to an eigenvalue integral. Further, for two sets of times we get again 
a normal matrix integral, where the eigenvalues of normal matrix fill the disc of unit radius \z\ < 1: 



P N (t,t;Y)=V^ 



db, 



UkiVk-bj) 



1 [dZ] exp V (ifcTr Z k + t k Tr Z f ' 

Ak,M<i Vfci V 



Tr ( B-\ A' +- (logZ+Z 



(14) 



These four boxed formulas constitute our main result. With the help of obtained expressions for the 
propagator one can construct matrix integral representations for different partition functions (|3|). While 
we do not advance too much in this direction, we present several examples related mostly to 2d YM 
theory: a "pants" amplitude, a genus one partition function and an expression for the simplest Wilson 
loopt 



2 Let us stress, that our partition functions does not literally coincide with the counterparts appeared in U(N) and 
SU(N) 2d YM theories. The ranges of the summation and the Casimirs are different. For example, SU(N) irreps are 



The structure of the paper is as follows: in Section Q] we remind the reader some basic facts about 
Schur functions and matrix integrals. In particular we remind the well-known orthogonality properties 
of the Schur functions with respect to integration over unitary and complex matrices, which allow one 
to glue different sums with each other. Then we construct three different representations of the Casimir 
operators. We derive the general expressions for operators of eigenvalue type and exponentiate them. In 
Section [5] we rewrite an exponential of the second Casimir operator as a matrix integral and act by this 
matrix integral valued operator on the initial conditions. In this way we get a propagator dependent on 
one set of times and on the matrix, which we marge with another function into propagator dependent on 
two sets of times. This propagator is naturally represented as normal matrix integral with the "square of 
logarithm" potential. Let us stress that the obtained in this section matrix integrals give the formal series 
representations of the sums over partitions. To obtain integral representations for convergent sums, which 
appear in some simple cases, one should consider an analytical continuation of the matrix integrals (jlOp 
and (jlip . In our paper we do not consider this type of continuation except for the particular example 
of the genus one partition function in 12.61 In Section [3] we repeat procedure of Section [2] for Miwa 
parametrization of the Casimirs coupling constants. In Section|4]we discuss obtained results and possible 
directions for the further investigations. 

1 Notations and basic formulas 

In this section we remind the reader some basic facts about Schur functions and matrix models. After 
that we introduce three different representations for Casimir operators, with one of them playing the 
crucial role in the subsequent derivations. We denote general matrices by bold capitals, for example X, 
while leave ordinary capitals for diagonal matrices. 

1.1 Schur functions 

The basic ingredient of the random partitions is general GL(oo) characters, which depend on infinite set 
of independent time variables tk and are labeled by a partition A: x\(t). Sometimes we put these times 
to be Miwa variables tk = jrTrX , and we freely denote this as a dependence on the matrix X, namely 
Xa(X) := xa(^/c(X)). For simplicity, all matrices are of the size N x N. 

Representations of GL(N) are parameterized by partitions A with the weight |A| = X^« an d the 
length l(X) <N: 

Ai > A 2 > . . . > Aj( A ) > - Ai (A)+1 = . . • (18) 

Explicit expressions for characters are given by Weyl formulas: either as determinant dependent on times 
tk 

X\(t) =detp Xi - i+j (t) (19) 

1,3 

where the Schur polynomials pk (t) are given by 

expl^ fe z fe ) =J2Pk(t)z k (20) 

\fe=l / fc=0 



labeled by the Young diagrams with Ajv = 0, so that 



>r W (U,V)= £ XA(U)x,(V)exp(^f^> +g |A|' 



(15) 



where 



cfW=£>(Ai-2i + l) + iV|A|-H_ (16) 

i — 1 

It is the term |A| 2 what breaks the Toda lattice integrability, which can be restored only in the limit N — > oo. For U(N) 
irreps are labeled by the "Young diagrams" without positivity restriction on the lengths of the lines, oo > Ai > A2 > . . . > 
Ajv > —00, so that annulus amplitude can be represented as a sum over representations with one additional variable r 

P^(U,V)= £ detlTdetV £ Xa (U) Xa (V) exp (| (c* u ™ + ^ + ^ ) + g|A|) (17) 



or, in Miwa parametrization, as ratio of two determinants 

Xa(X) 

where A(x) is the Vandermond determinant 



det x X ' +N ^ 

uel ',3 x i (21) 

A(x) 



A(x) = Y[(x t -xj)= dot xf 3 ( 22 ) 

i<j 

The very important role is played by the dimensions of the labeled by partition representations, of the 
symmetric group 

dim\(S N ) -i-r Xi-Xj+j-i 

dx = xx(tk = 5 k ,t) = W] = 11 f— 

0<i<j<oo J 

= TT (N- k)\ yr Xj - Xj + j - i 

J-J- (iV + A fc -fc)! H j-i 

and the general linear group 

dim A = XA (l) = dim A (GL(AO)= II ~ ~ ^ ~ ^' (24) 

0<i<j<A r ^ 

Here 1 is an identity matrix. The difference between dim A and d\ is that the first one explicitly depends 
on N, while the second does not. This difference is not so important for infinitely large N, but may be 
significant for finite N. The ratio of two functions is 



dim 



=n {x \t N _7v )l = NW ( i+ °i N - 1 » 



h i\ ( N -*y- 



(25) 



Let us also remind here the Cauchy-Littlewood identity 

exp ( JTfctfctfc ] = 5^xa(*)Xa(*) (26) 



\fe=l / A 

which, for the Miwa parametrization, leads to an identity 



exp ^i fc TrX fe = ^Xa(*)Xa(X) =: P N (t,X;0) (27) 



\/c=l 

and 



nl JV(JV-l) 1 
— —- = (-!)-»- det- 



1,3 



Xi-yj ' A(x)A(y) i,j Xi—yj 



(28) 



1.2 Matrix integrals 



In this paper we use matrix integrals with different integration domains and measures, in particular 
integrals over ensembles of unitary, Hcrmitian, complex and normal matrices. We denote these ensembles 
respectively it, Sj, <t and *Tt. In this section we fix our notations for matrix integrals and remind the 
expansion of the standard matrix integrals into the sums over partitions (character expansion). In all 
listed examples no Casimirs appear in the sums. Formulas of this section are scattered over different 
texts on symmetric polynomials and matrix models, for example J21 1l32f[34] , 

Unitary matrix integral. The most basic and important for our purposes is an integral over unitary 
matrices. We use the Haar measure such that the integral over unitary group is equal to identity: 

/ [dU] = 1 (29) 

JiX 



Then the following integration rules for characters are well-known 

/ [dU] XA(UAUtB) = Xa( ^ )xa(B) (30) 

JU dlm A 

J [dV] XA(UA) XM (UtB) = ^j^W (31) 

With the help of the Cauchy-Littlewood identity (|27p this leads us to the well-known matrix integral 
expressions for the propagator without Casimirs P]y(t,t;0) and for one of the vertices (J5]): 

L 

/■[dUlexpfg^TrCUAUtByU £ X*(*)xa(A)xa(B) ^ 

which are known as unitary matrix model and (generalized) Itzykson-Zuber integral. Original Itzykson- 
Zuber (IZ) integral for diagonal matrices A and B is a simple combination of their eigenvalues: 



[dU] exp J2 ^ Tr U fc + i fc Tr Ut fc = £ Xa(*)xa(*) (32) 

\fc=0 / 1{X)<N 



J [dU]expTr (UAU f B) = [ ][J fc! J 



detp aibj 

(34) 



A(a)A(6) 



In what follows we will mostly work with eigenvalue integrals. For example, an orthogonality condition 
([3i"j) in terms of eigenvalues for trivial A = B = 1 reduces to: 



N 1 r d 
Ilrf — |A(w)|2 XA(«)X^(S) = ^!^, M (35) 



i=i 



l«il=i "J 



Complex matrix integral. For the complex matrices we use a standard flat measure normalised 
by the constraint 

/[dZ]e- TrZzt =1 (36) 

and similar to (|3T?|) and (|3"Tj) integration rules for the family of the complex matrices are as follows: 

[dZ] e- TrZzt XA (ZAZtB) = XA(A 1 XA(B) (37) 

d\ 

f [dZ] e- TrZzt XA (ZA)x M (ZtB) = 2^5)*a,„ (38) 

A complex matrix can be decomposed in to the product of Hermitian H and unitary W matrices: 

Z = WH (39) 

where Hermitian matrix can be further diagonalized 

Z = WMU* (40) 

with unitary U and W and a real diagonal Mo Then the measure can be factorizedQ 

N 

[dZ] = N\v 2 N [dU] [dW] A 2 (0 n dl % (41) 



3 Let us note that real elements of M do not coincide with the eigenvalues of matrix Z, which are complex numbers. 

4 Accurate counting shows the discrepancy between the numbers of the degrees of freedom in the r.h.s. and l.h.s of this 
equality, namely 27V 2 real variables for complex matrix and TV + 27V 2 for Vs and unitary matrices. This is due to absence 
of U(1) N from the Cartan subgroup in one of the unitary matrices. As usual in the texts on matrix models division by this 
subgroup is assumed when necessary. 



where k are the eigenvalues of the Hcrmitian matrix ZZt = WAf 2 W^. We have introduced the notation 

N 1 

n * ( 42 ) 



V N 
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for the constant that is proportional to the volume of U(N) and is widely known in the theory of 
matrix models (see e.g. 35 ). Again, with the help of the Cauchy-Littlewood identity $Zjo§ the following 
expansions 



I [dZ] exp ( -Tr ZZ f + Y^ t k Tr (ZAZ t B) fe j = V 
Jc \ fc=i / i<\)< 



Xx(t)xx(A) X x(B) 



(\)<N X 



f [dZ] exp ( -Tr ZZ* + JT t k Ti Z k + t k Tr Z^ j = ^ 

^ V fc=l / ((A)<. 



XA(*)XA(t)dim A 

d A l j 

(A)<AT 

can be derived. Hermitian matrix integral. For the Hermitian matrix model we fix the integration 
measure as follows 

J [d9]exp(-1i^\=l (45) 

A Hermitian matrix can be decomposed into the product 4> = UJU* with unitary U and real diagonal 
X. The element of the volume is as follows: 

N 

{d<*>} = -^ K [dU}A(x) 2 l[dx l (46) 

( 27r ) 2 i=i 

It is simple to find a character expansion of Hcrmitian matrix integral. Namely, let us expand both sides 
of the identity 



/ [d*] exp ( -Tr \- Tr *Y J = exp ( Tr 



(47) 



in Schur functions of the matrix variable Y. With the help of IZ integral this gives 



$2\ XA^^Jdim* 

and, finally 



^ d*] exp (-Tr — j XA (*) = V ^ ( 48 ) 



[d*]exp -Tr— + 2^t fe Tr*' t =2^ -^ (49) 



V 2 k=l " " ' » dx 



The same expansion can be derived from the expansion of the complex matrix model (|44[) . Further 
we will use a matrix-valued delta-function 

/ [d*] exp(iTr *H) = 5(H) (50) 

which main property is 

/ [d*] *(# - H)/(*) = /(H) (51) 

for all (not necessary U(N) invariant) functions of the matrix variable /(4?). 

Normal matrix integral. This time integral is over normal matrices that is over matrices commu- 
tating with their conjugate, [Z, Z^] = 1. As usual, we fix the norm 

f [dZ]exp(-TrZZt) = 1 (52) 



A normal matrix can be diagonalized 

Z = UZU' (53) 

with the help of the unitary matrix U and the diagonal matrix Z with complex entries. Then the measure 
is 

N 

[dZ] = C £ [dU] |A(z)| 2 Y[d 2 Zl (54) 

For the normal matrix model expansions in Schur functions are tightly connected with those of complex 
matrix integrals, in particular 

/[dZ]exp(-TrZZt + f>Tr(ZZt)*)= £ ^hM (55) 

J<n fe=l 1(\)<N A 

Finally, we see that for constructing a matrix integral representation of the general sum ([3]) without 
Casimirs it is enough to use the vertices (1551) . (|4*5)l and the orthogonality condition (|3"Tj) . 

1.3 Casimir operators 

As we have seen in the previous subsection, classical matrix integrals give the possibility to construct 
matrix integral representations of the general sums over partitions Q without Casimirs. Thus the actual 
problem is to insert Casimirs into the sums. Casimirs (shifted symmetric sums), which we use 

do not coincide with actual Casimirs of GL(N),U(N) or SU(N), see e.g. [35]. The algebra of cut-and-join 
operators (Kerov algebra) is generated by Casimirs (1561) . Let us introduce operators 

C kX x = C k xx (57) 

There are at least three different representations of the operators Ck '■ in terms of derivatives with respect 
to time variables tk, matrix X or matrix eigenvalues Xi. For example, for the first (C\ = |A|) and the 
second {Ci = ^2 ^»(^» ~ 2i + 1)) Casimirs we have the following expressions in terms of eigenvalues^ 

A d 



^ = E^^T 



dxi 



N 



full matrix 



q = ST- X 2^_ y^ XjXj I d d_ 

2—1 ' i^j j \ j 



d 
Ci = TrX, 



(58) 



C 2 = Tr j X-— - ! - - ATi'X^s = TV : ( X 



d \ „„ „ d „ /„ d 



and times 



dX T J dX T ' V dX T 

8 



(59) 



Ci = j2 kt * 



fc=i 



dt k 



00 fi 
C 2 = 2J kmt k t m — h (it + m)t fe+ 



(60) 



dtk+ m dtkdtn 



5 Here we use the same notation for the operators of all three types, because further on we use only the operators acting 
on the eigenvalues. 
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Most convenient for our purposes are the operators that act on the eigenvalues. Let us prove a simple 
formula for the general Casimir operatorsQ From (|2Tj) the general expression for the Casimir operators 
in terms of eigenvalue derivatives immediately follows: 



iV 



e-W&i**) A(x) -^ 



where 



By definition 



A(x) = 



A(x) 



defc-X"*-* 



v=±(-^y=(-^ 



(63) 



(64) 



(65) 



is a constant, for example 






C 9 ° 



N 3 N 

~y~ 12 



(66) 



Let us prove explicitly that characters (j2"Tj) are eigenfunctions of operators (|63|) with eigenvalues C k ■ 



N 



C kXx (X)-C° kX x(X)=A-'\x)Y, 



dxi 



dctX"-* 



A ! 



a-'wE 



i=l 



' dxi 



E(-D H n 






£(Vi + J 



(67) 



Xa(X) 



which proves the statement. 

The crustal property of the operators (|63]l is that they can be easily exponentiated. Let us denote 
Xi — e Vi , then 



D(s) = exp ^2 s kCk = exp I - ^ s fc Cj? 

fc=l \ fc=l y 



1 



' °° N Q k 



^r{k Sk k^) A{eV) 



(68) 



2 Second Casimir 

In this section we consider only the first and the second Casimirs. Second Casimir is not only the most 
important for applications, but also is the simplest nontrivial one. It corresponds to the differential 
operator of second order, that makes it rather simple to operate with. To single out variables s\ and 



Let us mention a generating function of all Casimir operators in terms of times t k 
where 



Z)tt**= g~ r (^# — : exp (g(*e J ) -■£(«)) : -1 



fc=i 



Ck x k 


1 


V 2-wi J z 


fc! 


ea — e2 




k{z)-- 


OO / 

= £(?« 

fc=l V ft 



i a \ 

z k dp k J 



(61) 



(62) 



A proof will be presented elsewhere. 
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S2 in this section we denote them by q and g/2 respectively. Further on we will freely omit the explicit 
dependence on them: 

l(X)<N 

As we have mentioned it is simple to restore the dependence on S\, thus, for simplicity we will drop it of 
all intermediate formulas and restore it only in final expressions. 

2.1 Propagator with one set of times 

Let us construct a matrix integral for the propagator -FV(£, X). The operator (|68[) simplifies to 



where we denote by Dq a prefactor 

(gN gN 

A(e^) 
The identity 



N 

(£<%). A, «p(|E^)^*) (TO) 



D „ . " P ^ i I (71) 



=- !?£r =n^£*-p(-^< + ^) < 72 » 



helps to convert the operator ([701) into the exponential of the first order operator 

exp (§&) = -Ar /" <^A (e* + ») exp]T LA - 1, 
V2 > (2ng)^J-oo f^\ dipt 2g 

The shift operator acts on the "bare partition function" (|27[1 in the simple way: 



Pjv(t,X) = exp (|(7 2 ) exp K]t fc TrX fe 

\fe=i 



Do ' ■ 



l N yA (e* + y) exp]T ff>e fe (^) - ±yf\ 



2ff 



P N (t, X) = V / d"j/A (e») exp V V i fc e 

A,exp(-f Ef-A^) 2 ) r°° f^ N N 

= { *,^n / ^ A (eB) exp I j2 + E W(3H) 



where 



(74) 



This integral can be simplified by the shift of the integration variables 

y % ->yi- ifii (75) 

namely 



(76) 



2/i \ °° 

W(y) = -^+i K --NJy + Y / tk.e kv (77) 

The eigenvalue integral above can be represented as a matrix integral. Here we restore the dependence 
on q via change of variables t k — > tke kq which is equivalent to shift of the integration variables yi. As 
usual for Generalized Kontsevich Model, IZ matrix integral (|M)l gives the following: 
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[Sk] 





Figure 3: The propagator PN(t, t; s) combined of Pjv(t, X; s) and Pjv(t, X; 0). 



P N (t, X) = ^^ L / d M (Y exp -Tr $Y + Tr W Y (78) 

Here a matrix $ is diagonal with the eigenvalues ifi, and the integral is over Hermitian matrices Y. The 
potential is a minor deformation of W: 



N - 1 y 2 : ( q N 

and the integration measure is as follows 



W(y) = Wfo) + —-V = - y - +11- \ y + J2t k e k y (79) 

- \9 J k=1 



[d/i(Y)] = A(y)A(e») [dU] ne-^"'^ = exp ( £ ivWM ^^^ ) MY] 



(80) 



sinh ( 



Y®1-1®Y1 



= W det /TO1-1OT^ I dY l 



where [dU] and [dY] are usual measures for the unitary and Hermitian matrices, described m il. 21 Double- 
trace potential with coefficients made of Bernoulli numbers appeared in the non-flat measure (|80[) is 
similar to one effectively generated in the decomposition formulas [37]. Actually, it would be generated 
even for the simplest example, i.e. for decomposition of Hermitian matrix model, if one would know a 
matrix model or a simple field theory representation of the form tk = ( ex P(Sfc tk&k)) f° r the Kontsevich 
tau-function. A representation of this type is still lacking. 

Let us make a simple check of the obtained result. For N — 1 the definition (|69|) gives a simple sum 



Pi(t,x) = J2Pk(t> h e qk+iHk - X) (81) 



fe=0 



while the integral ([75)1 gives 



y 2 ( q + log x 1 



Jty 

fc=i 

3. OO ,m / 2 



e^ b 2 S / dyexp -— + I i 1/ -> /;,.-■ 



dyexp - + o * + 5> 



^ fc=1 



OO /'OO / 9 / 1 \ \ 



(82) 



2.2 Propagator with two sets of times 

The partition function Pjy(t,t) can be obtained from P/v(£,X) and a "bare propagator" expJ^tfcTrX 
with the help of the unitary matrix integral (|32p : 



p OO 

P N (t,t)= / [dVlPv^V^expVi/cTrV* (83) 

■to fc=i 
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Here we use an eigenvalue version (|35[) of the unitary integral. Combining it with (|74p one gets 

1 N 1 r rl / oo N \ 

W>^II^-f -^|A(.)| 2 exp EE^ & MM*.*) (M) 

where X is a diagonal matrix with eigenvalues complex conjugate to X{. One can omit the constraint 
\xi\ = 1 and integrate over a circle of the arbitrary radius Ri with all Xi substituted by x~ : 



x exp 



P N (t,f) = cTT -L / ^if d w t/A(:z)A(:z-V)x 

E (^ 2 + (jv - 1) to) - E (w^ ] + '**?))) 



(85) 



where 



eX P I 24 6 



(86) 
(27r 3 )tiV! 

Now we interchange the order of the integration with respect to a; and y and put Ri = exp ^ . Then one 
has 

and 



a^ 1 exp(j/ i ) = ajj (87) 

-oo - =1 ■ i "' J|xi|=e*3" '''J 

(88) 

E (^CiogN 2 ) 2 +(^- ^) togN a -E (**«? + **»?))) 

Z — X \ K/ — X / / 



Let us consider the integrals in Xj and y^ as integral over a complex plane C 



dmf ^~ i _^i (89) 



t k \ 



(90) 



which gives a normal matrix model (here we restored q via reorganization of Z): 

P N (t, t) = V- 1 f [dZ [, q exp ( --^Tr log 2 ZZ* + V (i fe Tr Z fc + t fc Tr Z 

where V is the normalization constant, which value is fixed by the constraint Pjv(0,0) = 1: 

V= f [dZ }e-^ Tll ^ 2zz "-( N+ i-^ l ^ zz ' =gf e £ N -% N2+ % N3 -* N f(N) (91) 

Jm 

This last matrix integral is both eigenvalue and Gaussian, thus, it can be easily evaluated. The key term 

in the potential is Tr log ZZ^; namely the terms of this form (while other details of matrix integral are 

different) are important in the matrix model for 3d Chern-Simons [1HB]. 

Let us make a simplest check of our result, again considering N = 1. In this case (|69|) is a simple sum 

of Schur polynomials: 

00 

P 1 (t,t) = E^(*K(*) e9 ' S+ffe(fe " 1) (92) 

fc=0 

while (|90l) gives 



. q 1 1 _£ 



Pi(^^- e " 2 y" ¥ /d 2 zexp(-llog 2 |z| 2 -Q-^log|z| 2 +^(^+t fc z fc )j^ 

(93) 

5 l_i_9 OO 



e 2 s 2 



a fc=0 



^ftWft(t) / di? e -*« 2 +( fc +f-i)« = una 

7 n J — OO 
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2.3 Specifications 

In this section we consider different specifications of matrix integrals (|78j) and (|§U)) and mention their 
relations to different applications. The connection with 2d YM we mention in this section is somewhat 
virtual: it would be precise for u GL(N)2dYM" gauge theory. For conventional SU(N) and U(N) gauge 
groups the sums and Casimirs can be (and actually are) different, see footnote [H 

1 . t — X propagator 

P N (t,X.) =J2xx(t)xx(X)e^ +S ¥ (94) 

A 

In addition to the Hermitian matrix model representation (j78[) this sum can be represented as a 
normal matrix integral with the help of the Miwa variables ik = -jjTrX fe . 

2. X — X propagator 

P N (X,Y) = J2 XA(X) XA (Y)e^ + i ^ 

l(X)<N 
r -2LTrlog 2 ZZ+-(Ar+I-.i)TrlogZZl- 

= V^ 1 / \dZ] - 

J m [ J det(l<g> l-X(g)Z)det(l(g)l- Y®Zt) 

This is a partition function on a cylinder for 2d YM and it can be used for the construction of 
higher genera partition functions as well as Wilson loops (see I2.6|) . 

3. t — 5 disc amplitude 

Pjv(Mfc,i)= E d *Xx(t)e^ + S ¥ 



(96) 



l(X)<N 

= V~ 1 f [dZ]e"^ Trlog2zzt " (A ' + ^t )TrlosZzt+TrZt+l: "=i*' cTrZfc 
Jm 

For infinitely large N this function gives a generating function of single Hurwitz numbers. 
4. t — 1 disc amplitude 

P N (t,l) = Y, dim xXx (t)e^ +S ^ 

l(X)<N 

r exp (-±Tr log 2 Zg - (y + \ - j) Tr logZZt + g^ t fc Tr Z*) 

At 1 J dfita-zt)^ 



jf [d/i(Y)]cxp [-^TrY 2 + («±I - ^ TrY + f>Tr 



e fcY 



5. Y — <5 disc amplitude 

sC 2 



P N (X,5 kA ) = Y,dxXx(X)e' lCl - 
x 

= -,-! /■ [dzl ^P (-^Tr log 2 ZZt - (Y + 1 - j )Tr logZZt + Trg) (gg) 

i m [ J det (1 (g) 1 - X ® Z) 



x / [d/i(Y)]exp( -Tr$Y-Tr 



Y 2 t fq N \ x; , , T , Y 



,Y + Tre 
/j, vy 2 5 V.9 2 

For q — this function can be considered as a generating function for the simple Hurwitz numbers 
in terms of the Miwa variables. Tn this case the last line of (j98| simplifies to 

66+21 /w., ^ / V-2A Ur 1 



¥ | d"j,A(* e ») exp (- £ | - (V - £ * + £ s ie ") (99) 



A(x)(2vr 5 )^ 
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This integral almost (up to substitution Vandermond A(xe y ) by A(y)) coincides with the matrix 
integral for the same function suggested in [27] . The identification of two integrals follows from the 
identity 



l -^ J d N yA( xe y) exp (- £ f g - (n - i) E Vi + E x t 
-|^/^A(,)e X p(- E |-(^-i) E , + E; 



(2^<r 

r\/A l-\i t m + vArln/itiAn /~\\ 4- Yi t~t nT\Avf»f/~iy /\ I 



(100) 



which can be easy proved by introduction of the operator A(^-) and consequent integration by 



parts. 
6. X — 1 disc amplitude 

_sC 2 



P N (X,1)= E dim AX A(X) e « Cl+i 

1(\)<N 

e -^Trlog 2 ZZt-(Ar+I-|)TrlogZZt " ' 



= 7>~ 1 / [dZ] . 

Jot det(l®l-X<g)Z)det(l-Zt) A ' 



This is the disc amplitude for 2dYM. While Hcrmitian integral (|78|l can also be simplified in two 
different ways, we omit here explicit expressions. 

7.6 — 5 spherical partition function 



P N (S k .uS k . 1 )= E d\e^ + S ¥ 

1(\)<N 

= V~ 1 f [dZlexpf - — Trlog 2 ZZ t - [N+- - - ) Tr logZZ* + TrZ + TrZ t 
Jm V 2 5 V 2 5 y 

For N — oo this should be equal to the partition function of CP 1 model [7H§] with only two first 
times switched on. 



8. 5—1 spherical partition function 



.Pjv(<$fc,ij 1) = E dimAdA< 



9C1 + 



.-/< 2 



, exp (-^Tr log 2 ZZt - (iV + 1 - |) Tr log ZZt + Tr z) 

At det(l-Zt)* 



/ [d/x(Y)] exp ( -^-Tr Y 2 + ( '±— - — ) TV Y + Tv , "" 



/c \ -1/ v g 2 



9. 1 — 1 spherical partition function 







Pjv(1. 


,l) = E dim 2 

A 


e lCi + ^ 










= p-i 


/ [dZ] 


e -iTrlog 2 ZZt-(AT+I-£; 

det (1 - Z) w det (1 - 


)Tr logZZ f 

-Zt) w 






e 


giV^giV.giV 2 

6 "■" 24 "T" 2 


q 2 N 

IT /" 


exp(- 


-fTrY 2 


+(*?- 


*) 


TrY") 




7V!( 27 r)f ff ^ 


A 


l«Ml * JJ 


det( 


1 - £*f 







(104) 



This is the partition function of 2YM on the sphere. 
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Figure 4: A genus one amplitude (lllip constructed of the propagator Pjv(Xt,X;g, g), 

2.4 Integrability 

Partition function (|90[) . as usual for matrix integrals, is a tau- function of the integrable hierarchy. To 
make a statement precise let us rewrite it as a determinant 

N 

P N (t,t;q,p) = V^W det Jnj(q,p) (105) 

where 

M«,p) = / tPzz*- 1 ^'- 1 exp f -1 log 2 |z| 2 - Q - ^ + i\A log \z\ 2 + £ (t fc z fc + i k J k ) J (106) 

It is obvious that hi y j(q + gN,g) does not depend on iV and 

dhij{q + gN,g) dhi,j(q + gN,g) , . 
^—-- =hi +k ,j{q + gN,g), *— -= = h iJ+k (q + gN,g) (107) 

This guarantees (see e.g. [35]) that the sum 

rjv(t,F)= ^ X AWxA(t)exp(^ (? (A l + 7V- i ) + |(A 4 + 7V- i ) 2 j (108) 

l(X)<N \i=l ' ) 

is a Toda lattice tau- function with respect to times £&, £& for arbitrary q and g, where N plays a role of 
the discrete time. 

2.5 ->■ limit 

Let us show that in the limit g — >• the partition function (|90p actually turns into the unitary matrix 
integral ([2]). Let us parameterize eigenvalues of matrix Z by real R and as follows: Zj = e~2~ +, "^> . Then 
in the limit 5 = one gets a product of the delta- functions ]X=i <K-^i — <?), which makes the integrals 
over Ri trivial. The resulting integral 



TV , / co N 



lim/\(/. /, = _[ j_./. __|A(,ir ; <-x P j > j >^ h e^+t fc e*tt5,) (ID!») 



is the eigenvalue representation of the unitary matrix integral (|32p . 

2.6 Synthesis 

Let us use the obtained matrix integrals to construct matrix integrals for some partition functions of the 
form (|3|). The first example is a three point function, matrix integral representation for which immediately 
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follows from the construction of the propagator ([TO]) : 

4"(M-Y) = Y XA(*)XA(f)XA(Y) e , Cl + fC2 = 

' ^— ' dim a 

/ A oo \ ( U °) 

= ^ / MZ i +1 i exp ( -1-IV log 2 Z Z t +f>Tr(YZ) fc +f fc TrZt*) ) 

•Art (detZZt)^ 5 * V 2 5 ^ v 'y 

Next important example is the genus one partition function, which is known to posses particular 
modularity properties [12) H5] 

z i°n(<1,9) = Yl ex P(9Ci + 2^2) (111) 

l(X)<N 

To obtain a formal matrix integral representation for this partition function one can just close the prop- 
agator (jnni) 

z { °$\ q ,g) = [ [dU]P w (u,ut) 

Ju 

r r exp (- -^Tr log 2 ZZ* - (iV + ± - |) Tr log ZZ*) 

= 7? Ju^ L [dZ] det (1 i 1 - U ® Z) det (1 ® 1 - IJt ® Zt) (U2) 

= 7>- 1 / dZexp j-^Trlog 2 ZZ t - (n + i - -J Tr logZZ 1 " + ]T ^TrZ fc TrZ tfc j 

In the eigenvalue representation of this integral all Vandermonds disappear due to the identity (|28[) . 
Obtained integral representation is suitable for investigation of the formal power series properties, but 
not for the investigation of the convergent series (jllip for negative g. To investigate this case it is more 
convenient to use representation (|74j) : 

N / oo 



(113) 



then 



( 



Z^(a a) - P V 2 6 24 ) 



N 
X 



n 1 I=I * £ * - (- g (I + ( w - 5) «* + ") ) as i^^ 

The last example here is the matrix model representation for the (generating function of) Wilson 
loops. For instance the simplest Wilson loop on the sphere is given by 

(Wx(C)) =^dim M dim, f [dU]x M (U) X ,(Ut) XA (Ut) 

w Ju (115) 

x exp (-^C 2 (/i) - ftdO*) - fc 2 (^) - gaCi(i/)) 
and the generating function of all such Wilson loops with generating parameters r, is as follows 

Z(r;g uqi ,g 2 ,q 2 ) = ]T XA(r)(W A (C)) 

i(A) ^ (116) 

- / [dU]Pjv(U,l;-si,-gi)P JV (U t ,l;-^;-g2)eS^ r '= TVUt " 
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3 All Casimirs 

In this section we turn on all Casimirs by introduction of Miwa variables for correspondent times 

s k = l T * Y ~ h ( 117 ) 

where the size of the matrix Y is M x M. For Miwa variables the propagator looks as follows: 

/oo \ ^ N M . _ i 

P N (t,i;Y)= J2 XA(i)XA(i)exp K]-TrY- fc C fc U ]T Xx(t) X X®UIl _\ ^ ■ 1 i 

;(A)<w Wi / ;(A)<w i=ij=i^ Al + * 2 

(118) 

3.1 Propagator with one set of times 

In this subsection we put N — M. We stress that this constraint is imposed just for simplicity of the 
description and can be easily omitted (for example, one can take some eigenvalues of the matrix Y to be 
infinite to effectively reduce its size). 

First of all, let us simplify the operator (p5|) , which for Miwa parametrization (|117[) looks like: 

N 1 

D(Y) = D J] __g-A(e*) (119) 



i, j= i y* d V] 



where 



D = A- 1 (enH(y j + i-l) (120) 

Using a complex matrix integral one can exponentiate the differential operator: 

D(Y) = D J [dZ] exp ( -Tr ( ZZ f Y - Z f Z^- J J A(e v ) (121) 

where •§- is a diagonal matrix with entries -£— . Let us act now by this operator on the "bare" partition 
function 

P N (t,X;Y) = D(Y)ex V (JTt k TrX k \ = 
= D / [dZ] exp ( -Tr ( ZZ f Y - Z f Z— J J A(e^)exp [ ^< fc TrX fe J = 

= D Q J [dZ]exp(-Tr (ZZtY)) f[ (f°° dA k 5(A k - (tfZ) kk )e Ak ^ A(e^)exp ff] t fc TrX fc j 

(122) 
After substitution integral representation of delta-functions S(x) — w- f™ e tpx dp and shift A k — > A k + ip k 
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one gets 

Piv(t,X;Y) 



^_ A [dZ]e -T r (zzt Y ) jj / /■ ^ / ^ufcC*™-^*)^ A (e A )exp ^, fcTr( 
= Ag^O jf [dZ] ^ [dU] JT [rfW] -q ^|~ ^ |~ ^ A{B) 2 A{A)A{e A )x 

x exp ( -Tr ZZ f Y + tTr UBU+ (WiW 1 - <p - Z f Z) + ( - - N J Tr A + ^ < fc Tr e* 

= A) AM / [dZ] / [dn(A)} f [dB] 

( N °° \ 

x exp -Tr ZZ f Y + iTr B (A - (p - Z f Z) - — Tr A + ^ t fc Tr e fcA 



\ fc=i 

where we introduce 



(123) 
A = WiW f , B = UBU 1 (124) 



and use the definition (|80[). Then one can integrate out matrices B and A using representation (|5U)l of 
the matrix-valued delta-function to get 

Pjv(t, X; Y) = D A(<p) [ [dZ] exp (-Tr ZZ^Y + #(Z+Z + <p)) (125) 

where the potential 

H(A) = -^TrA + f>Tr e - + £ ^fff Tr A * Tl ' A ' (™) 

fe=l i;j"=0,i+j>0 V J ' '■'' 

Let us stress that since we have integrated in angular degrees of freedom in the nontrivial way, the 
obtained matrix integral does not directly simplifies to a eigenvalue one. Let us check the consistency of 
the derived matrix model, namely to check it for N = 1 and, perturbatively in t, for arbitrary N with 
X = l. 

For N = 1 the propagator (|118[) simplifies to 

oo 1 

Pi(t,e*;g)=X>fc(*) e * y Hi (127) 

fc=o y fc + 2 

In this case 

C23= (?/+i) e¥ / dmexpl-my-^(m + ip)+J2 t 3 eKm+V) )^C22il (128) 

For X = 1 (ipi — 0) the integral ( |125|) simplifyes to 

N 



P N (t,l;Y)= f[ (iVj+^-^x 
JjdZ}exp (-TrZZt (y + y) + f] i fe Tre fcztz + £ i=i^^Tr (^Tr (ZtZ)' 



(129) 



20 



Using an explicit expression pTj) for the measure on the space of complex matrices one easily integrates 
out angular variables with the result 



Pjv(t,i;Y) 



el,=i fa + * 



- = ^ A(-y) / d^ACe^dete-^+^-^^exp VVt/' = 

I) n*=i ^o « v^ 1 ^ / 

WAr7V! r^dete-^ +i -^exp(f;f;t fc e^] = (130) 



A(-y) 



/ i N 

det r + ti V det 

\* -i =1 1A + J - 3 ^ 4 J' =1 



It is easy to see that in the last sum all terms except for the first one are equal to zero. Thus, with the 
help of (|28|). one gets 

N + - 
P N (t,l,Y) = l + t 1 N]J^—f + ... (131) 

3=1 V] 2 

which coincides with the first terms of the expansion of (I118[) for X = 1 . 

3.2 Propagator with two sets of times 

To derive the propagator with two sets of times one can use the propagator with one set derived above, 
but here we use a slightly different approach. Namely, we apply the following representation 

nn-V=niW *** ** , h] ^ 

where a contour C encloses all poles of the denominator of r.h.s (thus the contour integral gives just a 
sum of residues in the points y\ . . . , tjm)- After substitution of this operator into (|119|) one gets 

N i r r°° e a i h i ~ ( °° N \ 

P N (t,X;Y) = D l[—idb ] / dOj-sr- —A(e*-°)e^2> $>*<<"-*> (133) 

j=x *mj Jo Uk=i(yk-bj) y k=1 j=1 J 

Then, to get the propagator dependent on two sets of times, one can apply the same trick as in the 
previous section ([55"j) -(j!TT j) with the result 

N 

P N (t,t;Y)=V Y - 1 d 2 z i |A(z)| 2 JJexp(^(z i ,zi)) (134) 

J\Zi\<l l = 1 

where the potential is given by 

exp (W(z, z)) = I db-j^- exp ( f] t k z k + t k z k - f^ (b + N + M log \z\ 2 I (135) 

J c n fc =i(2/fc - fo ) \ k =i j=i v l) ) 

In terms of integrals over matrices this transforms to 

P N (t,i;Y) =Vy- 1 <fdb J TJ , l - f [dz]e (l2^ k ^ + t k T^)-T,(B + N +i ) los ^z) 

JC [[kiVk-Oj) Jffl,\zi\<l 

(136) 
Of course, this is just a representative of the possible matrix models, for example, one can absolutely 
similarly construct a normal matrix integral with eigenvalues situated not inside but outside of the circle 
|z| = l. 
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One can integrate out all bk to get 

N M / -, \ M M 



iV IV1 / -. \ 1V1 IVJ -. .. 

p * Y »~nn » + <- E-i n ^-- n ^ 

1=1 J = l v ' Ji=l ljv = lfci5tli 1 A k N =£l N " " 

„ / oo JV / ^ 

x / [dZ] exp V (i fc Tr Z fc + t,Tr Zt fc ) - V ( w + TV + - 



x 



(137) 
loe;|z,- 



Let us make the simplest check of the obtained matrix model representation. For N = 1 and arbitrary 
M expression (I118[) simplifies to 

oo M i 

Pi(t,f,») = £p*(*)p fc (t) n ym ,ii ( i3s ) 

fc=l m=l ^ m ~ fc + 2 

Formula (|137[) for N — 1 gives: 

A/ / -, \ m / \ „ / , „ x oo 



+ > fez fc + £ fc z fe ) 



n(«^)E(n i ^JX 1< /-p(-(« + I)v<i.p + E(« 

;=i v y i=i y m ^i ym y '/ fe=i j i^i< i 



(139) 



We do not discuss different specifications of the obtained matrix models here. Let us only mention 
that the matrix model (|137l) for tk — ifc = <5fc,i is similar to one conjectured for (stationary sector of) 
the CP model in |18) , but our model is more involved. Expressions for different higher genera partition 
functions, in particular, for genus one, can be constructed as in previous section. 

3.3 Integrability 

Again, as in the previous section, we present the propagator (|136j) as a determinant 

N 

P N (t,t;Y)~ det h UJ (140) 



where 



/ db^— I Szz 1 - 1 ^- 1 exp [ V (t 

Tc U k (yk-b)J lZil<1 \t^ K 



hi,j = j> db^^ ^ / d 2 ** 4 - 1 ^- 1 exp [ Yl (t k z k +i k zJ ) - ( b+ - ) \og\z\ 2 ) (141) 



tfe\ / , 1 

2 

"«-*— - - 1-1^- \k=\ 

Equations 

— = hi- k!J , — = hi!J _ k (142) 

guarantee Toda lattice integrability of the sums 

CN 00 \ 

^2J2 Sk ^ +N -^ k ) ( l43 ) 

with respect to times i and i. For N — > 00 the sum in the potential can be regularized as in [TU] 
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4 Conclusion 

In this paper we construct the precise relations between random partitions of the finite size and matrix 
models. A huge number of interesting topics, such as the phase transitions of the obtained matrix 
integrals [21] , their role in M-theory of matrix models (decomposition formulas) [37j and Virasoro-type 
constraints [23[3S] (in particular, an application of the powerful Eynard technique 22,25,26 ) for them 
are beyond the scope of this letter. 

It is not obvious to us if there exists an infinite set of the Virasoro-type constraints that is an algebra 
of the low-order differential operators, which act in the space of the time variables tk, tk and, probably, 
Sfe, and annihilate the propagator JV(i, i; s). At least for the matrix model (|90[) a usual invariance of the 
matrix integral does not lead to the constraints representable in the differential operator form, except for 
the eigenvalue rescaling Zi — > e e Zi, Zi — ¥ e e Zi, which leads to the equations: 

00 f) °° ft ft 

£^^ = £fe^ = ^ (144) 

fe=l K fc=l H 

These equations are obvious from the definition of the function Pjv as a sum over partitions. We guess 
that to obtain the closed set of Virasoro-type constraints one should be able to introduce additional 
observables into the model. 

As far as concerns integr ability: partition function ([7]) for infinite N is known [36] to be the Toda- 
lattice tau-function with times t, t. Both for the second Casimir (|90|) and for all Casimirs coupled with 
Miwa variables (I136[) this integrability is obvious from the matrix model representation for arbitrary 
finite N; and N plays the role of discrete Toda time (this can be also derived directly for the sums over 
partitions from the considerations in 36 ). Integrability of another type, namely Toda-chain integrability 
in times s& for tk — iu — $k,i EZnSI IT] - is by no means obvious from our matrix model representations. 

Operators Ck , which let us to construct non-trivial partition functions with the help of the basic ones 
are similar to the operators appearing in an under-developed theory of the check-operators [39] . 

Let us also mention one of the possible generalizations, which is extremely interesting, namely /3 
generalization, important the recent AGT conjecture [T9j[20] as well as for other applications. The 
problem with this generalization is that while we know very well a proper /3-generalization of Schur 
polynomials, namely Jack polynomials [291132] . /3 ^ 1 analogs of Casimirs operators are not so simple to 
operate with. For example, the analog of the second Casimir is the Calogero-Sutherland Hamiltonian 

A-t«J5T'i;^(s;-£) (145) 



and we did not manage to find for its exponential any simple analog of (1701) . 

In the subsequent publications we are going to consider a very interesting and important for applica- 
tions question of large N expansion and topological expansion. Let us just mention here, that in different 
applications of the random partitions appear two different types of genus expansion: 

1. In "2d YM"-like sums, when summands are combinations of xa(X), dim \ and Casimirs, the role of 
the topological expansion parameter is usually played by the natural for matrix models 4, constant. 

2. In "Hurwitz-Hodge-Gromow-Witten" partition functions, with sums built of xa(£) and d\, one 
usually does not consider A corrections and simply puts N — oo, that is the summation is over all 
representations of GL{oa) without any restrictions. In this case the topological expansion goes in 
additional parameter h, which we do not introduce in this note, see e.g. [?TUlj. 

In the last case the situation is in some sense intermediate between Hermitian matrix model, for which -^ 
plays the role of the topological expansion parameter and Kontsevich model, for which N just counts the 
number of the independent time variables and does not explicitly show itself in the partition function. 
We conjecture that ■h corrections even for "Hurwitz-Hodge-Gromow-Witten" partition functions contain 
important physical information and should be investigated. Here the simplest example is the unitary 



23 



matrix model 



Zg$ } (t, t; 0) = [ [dU] exp ( £ £ fe Tr U fe + i fc Tr ut fc j = £ XA ( <)xa (j) 

^ U \fc=0 / 2(A)<JV 



(146) 



which for any finite iV is a highly nontrivial function of times t, t, but for N = oo it transforms into 

4 ;2)(m;°)=£xa(*)xa(*) ( i47) 

A 

which, due to the Cauchy-Littlcwood identity (J26I) is a trivial exponential exp (y^,_ 1 ktktk)- 
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